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' We consider static spherically symmetric stellar configurations in Palatini the- 

ories of gravity in which the Lagrangian is an unspecified function of the form 
f(R,R fJiU R flu ). We obtain the Tolman-Oppenheimer-Volkov equations correspond- 
ing to this class of theories and show that they recover those of f(R) theories and 
General Relativity in the appropriate limits. We show that the exterior vacuum 
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solutions are of Schwarzschild-de Sitter type and comment on the possible expected 
modifications, as compared to GR, of the interior solutions. 

PACS numbers: 04.50.Kd, 04.70.Bw, 04.40.Nr 



INTRODUCTION 



One of the most challenging open questions that Physics faces nowadays is that of 
explaining why the Universe has the structures that we observe and how they came into 
existence. The current establishment states that some yet-to-be-determined source of dark 
energy should be responsible for the observed speed up in the cosmic expansion rate 
^ " Additionally, it assumes the existence of huge amounts of dark matter, which must have 

been necessary to increase the rate of growth of structures out of the highly homogeneous 
initial configuration observed in the cosmic microwave background radiation [2j (see also 
[3|). The effects of the dark matter component should be observable today in the dynamical 
and kinematical properties of stars in galaxies and of galaxies in clusters. 
Though the dark matter and dark energy models are able to successfully fit different data 
sets, they also suffer from some limitations that motivate the search for alternative expla- 
nations (see for instance In this sense, scenarios in which the gravitational interaction 
is modified at large scales have been proposed to address the same phenomenology as 
the dark matter and dark energy models. This has led in the last years to investigate 
theoretically and observationally the effects of possible deviations of the gravitational dy- 
namics in the Newtonian and Einsteinian regimes (see, for instance, the review articles [JB])- 
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Though the existence of dark energy has very little effect, if any, on the structure 
of stars, the presence of dark matter or of modified dynamics can have a much more 
dramatic impact which, in some cases, could lead to observational effects. In fact, it has 
been suggested that dark matter could play a non-trivial role in the mechanisms of energy 
generation and transport in stars, which could be used to set bounds on the type and 
properties (decay rates) of the hypothetical particles making up the dark sources |7| • The 
potential effects of modified gravitational dynamics are more difficult to estimate, however, 
because one must derive the precise set of equations that govern the internal structure of 
fluids in equilibrium and then consider specific models of interest, such as relativistic stars 
or compact objects, which unavoidably requires the use of numerical methods. Despite 
these difficulties, the study of the effects of alternative equations on the structure of 
compact objects, stellar formation and evolution, and peculiar objects (instability strips, 
protostars, etc) has attracted some attention in the context of modified theories of gravity 
of different types [8-11]. Given the increasingly high accuracy and resolution of currently 



available observational techniques (see for instance [12J, [13|) and the discovery of rare 
objects which defy the standard model of stellar formation and evolution [jij . the study 
of stellar properties and new solutions in extensions of GR could provide new insights to 
interpret observations and offer new avenues to test the strong field regime of gravitation 
and, therefore, help determine or set limits to potential corrections to Einstein's equations. 



Motivated by all these reasons, in this work we study how the structure equations of 
perfect fluids in hydrostatic equilibrium, the so-called Tolman-Oppenheimer-Volkoff (TOV) 
equations, are modified with respect to those found in General Relativity (GR) when the 
gravitational Lagrangian is extended to become a function of the Ricci scalar and the 
Ricci squared scalar, i.e., f(R,Q), with Q = g^g^R^Rap = R^R^. Though the usual 
metric formulation 1 of this type of theories is known to have ghosts and other dynamical 
perturbative instabilities due to the higher-order character of the field equations, we 
consider its Palatini (or metric-affme or first-order) formulation jBf, which is free from those 
drawbacks, as will be explained in detail below. The family of Palatini f(R,Q) theories 
have been studied in the recent literature specially in scenarios involving high-energy 
physics, such as early-time cosmology 15], black holes jl6j |. and applications to quantum 

WM. In 



gravity phenomenology [17], and also in the context of the cosmic speedup 
the usual (metric) formulation theories of that kind are also well known (see, for instance, 



2lL |22| for some early references, and |6j for more recent works). 



The Palatini formulation of f(R, Q) theories leads to a set of second-order equations for 
the metric that exactly recover GR in vacuum. In regions containing sources, such as in the 
interior of stars or fluids in general, the dynamics is modified due to nonlinear matter terms 
induced by the form of the Lagrangian. These nonlinearities arise due to the nontrivial 
role played by the matter in the determination of the affine connection, which is assumed 
to be independent of the metric (Palatini formalism), i.e., it is not constrained a priori 
to be given by the Christoffel symbols of the metric. As a result, material systems whose 
gravitational dynamics is governed by a Lagrangian different from that of Hilbert-Einstein, 
f(R, Q) = R, could admit equilibrium configurations different from those found in GR, 



The metric formulation assumes that the connection is constrained a priori to be given by the Christoffel 
symbols of the metric. 
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which could have observable consequences. 

Although Palatini f(R,Q) theories do not introduce higher-order differential equations, 
they involve elaborated algebraic manipulations (required to solve for the connection) that 
lead to nontrivial modifications of the dynamics. For this reason, in this work we focus on 
the derivation of the corresponding TOV equations and on putting them in a suitable form 
that facilitates their use in numerical applications, which will be considered elsewhere. We 
also compare the resulting expressions with those obtained in the literature for Palatini 
f(R) theories. 

The content of this paper is organized as follows. In Section [Til we review the formulation 
of f(R, Q) theories in the Palatini formalism. We then derive explicit expressions for the case 
of perfect fluids in subsection III Al and for simple models of the f(R, Q) action in subsection 
III Bl The stellar structure equations in those models are derived in Section III II where we 
also check that GR is recovered when the appropriate limit is taken. These results and its 
potential applications are discussed in Section IIVI 



II. ACTION AND FIELD EQUATIONS 

We define Palatini f(R, Q) theories as follows 

%,r,vg = 2^/ d A x^gf(R,Q) + S m [g^ m ], (2.1) 

where n 2 = SnG, S m [g,ip m ] represents the matter action, g a/3 is the space-time metric, R = 
g^Rfiu, Q = g m g v ^R i ivR a p, R^v = R p fJ , P u, and R 01 ^^ = d^T^ - d u T^ + r^r^ - T" X T^ 
Variation of (12. ip with respect to metric and connection leads to the following equations 23 



f a 2 

IrR^u — -^g^u + 2fQRnaR a u = « T^ v (2.2) 

V/3 [s/=g {Jsgr + 2/qJO] = . ( 2 -3) 

were we have used the short-hand notation fx = dxf- For simplicity, in the above derivation 
we have assumed a symmetric Ricci tensor, R^ = R„,, ., and vanishing torsion (for more 
details on the effects of relaxing these conditions, see [24|). As shown in [23j, the connection 
equation (12. 3p can be solved introducing a rank-two tensor (or auxiliary metric) h^ u such that 
V-9(fR9'" / + 2 fQ Rflu ) = V^hh^, which turns Q into [\^h^ v ] = and implies 
that the connection can be expressed as the Levi-Civita connection of h^ v . Therefore, the 
Ricci tensor of the independent connection is equivalent to the Ricci tensor of the metric 
h^, which from now on we denote R^ u (h). With simple algebraic manipulations, one finds 
that the relation between h^ u and g^ v can be expressed as follows 

V = V^ditEfE- 1 ]/^ , hT = y -j=^ , (2.4) 

Vdet E 

where we have defined the matrix T> a u = (fn^ + 2fqB a v ) and B a v = R a fi{K)g^ v . It is 
important to note (see section III Bl below for details) that the matrix B a u is a function of 
the components of the stress-energy tensor. This implies that R = B a a and Q = B a v B v a 
are functions of the matter. According to this, the physical metric g^ v and the auxiliary 
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metric are related by a matter-dependent deformation given by the matrix Tj a u . This 
deformation becomes a conformal factor in the particular case of Palatini f(R) theories 
(where /q — > 0), as is well known (5)]. 

With these relations and definitions, the field equations (12. 2p for the metric can be 
written in compact form as B fl a T la u = ^8 U + k 2 ^", and using the relation B/S/ = 
V / detS R lia (h)h aL ' we obtain 

v( ^7ab(^ + KV )- (2 - 5) 

This form of the metric field equations will be used in what follows to derive the TOV 
stellar structure equations. However, before that we need to determine the explicit form of 
the matrix Eq,^. That task is carried out in the next subsection. 



A. f(R,Q) theories with a perfect fluid 

The explicit form of the matrix £ that relates the metrics h^ v and g^ u can only be found 
once all the sources that make up T^ u have been specified. In our discussion we will just 
consider a perfect fluid or a sum of non-interacting perfect fluids such that 

T„ v = (p + P)u^u u + Pg^ (2.6) 

with p = Y2iPi an d P — TltiR- 111 order to find an expression for E, we first rewrite (12. 2p 
using matrix notation as 

2f Q B 2 + f R B - f -I = k 2 T . (2.7) 
Using ( 12. 6 p this equation can be rewritten as follows 

/•B A ( ,2 r> 1 / 1 ffl 



2fQ \ B+ W) = v + ip+ p)uX ■ (2 - 8) 

Denoting A 2 = {^k 2 P + f + jf^\ an d making explicit the matrix representation, (12 .8p be- 
comes 

fa A 2 (\ 2 -k?{p+P) 



where /3X3 denotes 3-dimensional identity matrix. Since the right-hand side of (I2.9P is a 
diagonal matrix, it is immediate to compute its square root, which leads to 

where si denotes a sign, which can be positive or negative, and S 3 X3 denotes a 3X3 diagonal 
matrix with elements {si = ±1}. For consistency of the theory in the limit /q — > 0, we must 
have s\ = 1 and S3X3 = -^3x3- This result allows to express the matrix E^" as follows 



ai 
02/3x3 



(2.11) 
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where a hat denotes matrix representation, and 0\ and a-i take the form 

°2 = y + • (2-12) 

Note that we have kept the two signs in front of the square root of a%. In order to correctly 
recover GR at low densities, one must take the positive sign in that equation. However, 
at high densities the square root may vanish and one may need to take the negative sign 
branch to guarantee that a± is continuous and differentiable accross the point where the 
square root vanishes (this subtlety in the behavior of o\ was first observed in [l5|). This 
technical issue does not arise for <r 2 . 



B. Workable models: f(R, Q) = f(R) + aQ 

So far we have made progress without specifying the form of the Lagrangian f(R,Q). 
However, in order to find the explicit dependence of R = B^ and Q = B^BJ 1 with the 
p and P of the fluids, we must choose a Lagrangian explicitly. Restricting the function 
f(R,Q) to the family f(R,Q) = f(R) + aQ, we will see that it is possible to find the 
generic dependence of Q with p and P, while R is found to depend only on the combination 
T = -p + 3P |T 



25 



The reason for this follows from the trace of (12.21) with g^, which 
for this family of Lagrangians gives the algebraic relation R/r — 2f = n 2 T and implies that 
R = R(T) (like in Palatini f(R) theories). For these theories, we have that /q = a, which 
is a constant. Therefore, from the trace of (\2.9\i we find 



which can be cast as 



X 2 -k\ P + P) 



After a bit of algebra we find that 



A 



/q 



/q 



From this expression and the definition of A 2 , we find 



(2.13) 



(2-14) 



(2.15) 



3 I R + ^ I + 



\ 



where R, f, and fn are functions of T = — p + 3P. 



r . h) 4^(p + P) 
/q 



(2.16) 
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III. STELLAR STRUCTURE EQUATIONS 



A. Geometric part. Preliminaries. 



Given the field equations in the form ( 12.51) and the matrix £ of a perfect fluid, see 
(12. lip , we have all the elements to compute the Ricci tensor R^ U (T) = R^^h), which 
represents the left-hand side of the field equations (12.51) . Choosing the diagonal metric g^ v 
as 9/j.v ~^ (— A(r)e 2 ^( r \ 1/A(r), r 2 , r 2 sin 2 9), the corresponding h^ u diagonal elements are 

h tt = J— 9tt = Sg tt (3.1) 
hij = y/(Tia 2 gij = Vtg i:j . (3.2) 

The inverse metric components are trivially found from these ones. The non-zero Christoffel 
symbols are (obtained by direct computation) 



L tr 
rr 

J- A.A. 



1 dr[Sgtt\ 

2 Sg tt \ 

: sin 2 6TI, 



pr 



pr 

rr 



dr[gttS] -p r 
1 8 

r 



2ttg 

dr [Qg r 



2ng 

For completeness, we expand those coefficients as follows 

v\ r = l(^ + 2A + f) r& = -llg + aw + f) r 



re 



2Ur 2 
cos 6 
sin 8 



r 2 A (Vt 



rr 



sin 2 01^ 
— sin 9 cos 



1 (Or _ Ar\ 

2 \ n A ) l t8 

rt = rL rt 



pr 

rr 



2 

1 (VL 
2 

cos 
sin 8 



(f + f) 



The terms that contribute to the Ricci tensor components are the following: 
R tt = d r T r tt + T r tt (T r rr + 2T e rg - T* t ) 

Rrr = 'dr (i* + 21*) + 1^ (i* + 21* ) - ((J*) 2 + 2(r^) 2 )) 

Ree = 1 + <9 r rg e + (Y rt + T r rr ) 
R^ = sin 2 9 R ee 

Inserting the corresponding Christoffel symbols, the result is 



Rtt 



A 

■n-rr 

T 



1 (Sg tt 

2 \Vg 

rr 

A r S r 



+ 



+ 



l 

2 



A 



rr 

~A 



A,- 
~A 



Re 



A ri S 
A(l + rip r ) - 



+ 2ip rr + 



A 



rA r 




(3.3) 
(3.4) 
(3.5) 
(3.6) 



(3.7) 



(3.8) 



7 



B. Metric field equations 



We now focus on the equations of motion (I2.5p . Using the notation R^ 
the following useful relations 



', we find 



O t 

—lit — ti*r 



R — T e 

tie — 

Using the results of above for the Ricci tensor, we find 



(3.10) 
(3.11) 



n 

s 



Rf — R r 



Rq 



-A 
2M, 



Ai\) r 



A 
~2 



1 1 2 r . ( _ 1 \ r St \ X ( Sr- - 2 r 

-I 9 - -I - 1 H I — - - 

20 \ H 5/ r o 

a, / 4 A r \ i /o r 2 



rr 



We can thus write the structure equations as follows: 



(3.12) 



2M r 



1 

re 6 



n 



2— + — 

2 o V 5 



.4 
2" 



O r /4 A 



i /a. 



o lr + i J + 2\ n + r 



1 ( Sj> 



+ 



+ 



o~ 



(3.14) 



The second equation can be further simplified using (13.141) and becomes 



K - §r< + A 



— 11 _| — H - 



(3.16) 



where we have used that for a perfect fluid rg 9 = t/ . With a bit more of extra effort, it 
can be put in its (almost) definitive form as follows 



2\ Mr 3tI 



O /■ / r 



s r t 



+ A 



Q rr Q r f 2r - 3M 3 O f . 
+ TT" Vr(r-2M) ~ 4~0 



(3.17) 



This last expression can be directly compared with equation (11) of [lO], where f(R) 
theories were considered. The f(R) limit is recovered taking Q = S = fn, and a wrong 
factor in 101 is corrected. 



C. Conservation equation 

Let us now focus on the conservation equation V^T^ = 0, which takes the form 
dP , s f A r \ , ■( M , M r \ 
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Combining f !3.14j) and (j3.17p . we find 



ft + r 



M r 
rA 



ft r T3ft r 2r-3Ml 1 fS r ft r \ /ft r 2 

2A "ft |_4 "ft + r(r - 2M)\ ~2\S~!T) \ tt + r 

(3.19) 

It is important to note that this expression is independent of ft rr , which guarantees that 
(13. 181) will become a second-order algebraic equation for P r , as it happens in the f(R) case 
and in GR. Denoting A = (^f + ~), Q5ZEB} becomes 



dP 

dr 



<P + P) 



M 



2AA 



1 ft r 
A "ft 



3ft r 
4 "ft 



2r - 3M \ 
r(r - 2M ) J 



1/5; 
2 1 5 



ft r 
"ft 



(3.20) 

Assuming an equation of state of the form p = p(P), we can formally denote ft r = ftpP r 
and SV = SpP r , where ftp = <9pft and Sp = dpS. With these definitions, (I3.20p can be 
written as 

d + C 2 P r + C 3 P r 2 = , 



where 



c 2 
c 3 



2(P + P) 
r 2 (r - 2M) 
2 



M 



S* 4 



r 
ft P 
"ft" 



(p + P) /ftp 

2 v ft s 



(p + P) (3 ftp /5 



2 ft 



5 



ftp 
"ft" 



(3.21) 

(3.22) 
(3.23) 
(3.24) 



One can now compare these expressions with those obtained in [iCj for the f(R) case. We 
find two typos, one of which corresponds to propagating the error previously found in (I3.17p . 
and affects the factor 2 in front of C%. The second typo appears in the definition of C3, where 
the factor 1 that appears in the square bracket is missing in [lol ]. Using the same notation 
as in 



10J, we can express the pressure as follows: 



>(0) 



[1 - a(r)] 



l±dl- p{r)P, 



(0) 



(3.25) 



where the + sign in front of the square root should be taken to recover the GR limit, and 
have used the shorthand notation 



p(0) 
r 

a(r) 
0(r) 



(P + P) 
r(r - 2M) 



r 

7 



(p + P) (Op Sp 

2 V ft s 



(2r) 



ftp 
"ft" 



(p + P) /3ftp 
2 ft 



ftp 
"ft" 



Sp 
5 



(3.26) 
(3.27) 
(3.28) 
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D. Summary of results 

We can now write together the three equations that determine the stellar structure in 
/ (R, R pil R liV ) Palatini theories of gravity: 



n r ) r 



1 

A 

3tI 



s T t 



+ A 



~~n + IT 



'1T + ^) ~r 
2r-3M 



r(r - 2M) 



3 a 

4~fi 



>(0) 



[1 - a(r)] 



l±Vl-/3(r)P, 



(0) 



(3.30) 
(3.31) 



E. Additional manipulations needed 

In order to put the equations in suitable form to do numerical calculations, we still need 
to specify the form of Q rr . Note that once an equation of state is given, we can express 
Q r and S r as QpP r and SpP r , respectively. The term Q rr needs some extra manipulations, 
because it has the form Q rr = QppP^ + QpP rr . After some algebra, one obtains the following 
expression for P rr 



1 rr 
~Pr 



P 



(0) 



p 



(0) 



1+ l 



(3P 



(0) 



(3P^\l + s\/l- (5P 



>(0)> 



+ 



a r s 
— — + - 
1-a 2 



(3 r P 



(0) 



where s = ±1, a r = apP r , f3 r = [3pP r , and [here we define $ = (r^ + §T f *)] 



i-^i + ^i-^; 

(3.32) 



,(o) 



P 



(0) 



l±pp 
p + p 



2(r — M) r s ( $ P P r + ^ 



Pr- 



r(r - 2M) 4 



r-2M M-$^ 

4 



(3.33) 



This last equation can be equivalently written as 



p(0) 



P 



(0) 



1 + Pp 



p + P J (M-$£) 



3$r 2 



2(r - M) 
r(r-2M) M - $2± 



+ 



1 



r-2M M-$^ 
(3.34) 



This shows that P rr can be expressed in terms of r, p(P), P, P r , M, and M r . Therefore, our 
system of equations only involves first derivatives of the functions P r , M r , and ip r . 



F. Limit to f(R) and GR 



Equations (I3.29p . (I3.30p . and (I3.3ip allow to obtain the TOV equations corresponding to 
Palatini f(R) theories by just taking Q = S —¥ fn and setting f(R,Q) — > f(R). One then 



10 



obtains 



/b r 



1p r 

M r 



1 



.4 

— £ 1 + A 



+ 



+ 



/s,r ( 2r - 3M 3 /p r 



/h \,r(r - 2M) 4 /, 



(3.35) 



The expression for the gradient of the pressure is the same as ( 13.311) except for the fact 
that Pr°\ oi(r), and /3(r) get redefined as follows 

~3" 



p(0) 



(P + P) 
r(r - 2M) 

fR,P 



M - ( r ; + rt) 



a(r) = (p + P) 
P(r) = (2r) 



Ir 
Ir,p 
Ir 



l-(p + P) 



3/ 



4/ 



(3.36) 
(3.37) 
(3.38) 



The limit to GR simply requires to take /r — >■ 1 and /(-R) = -R after taking the f(R) limit, 
which leads to 



r 

2M r 



1 

A 



/t 2 (p + P) 



.4 



3r: - r* 



2 

(P + P) 



K 2 p 



r(r - 2M) 

(P + P) 
r(r - 2M) 



M-tf + rt)- 



M + 



K 2 Pr 3 



(3.39) 
(3.40) 

(3.41) 



where A = 1 - 2M(r)/r, and in GR t/ = n 2 (p - P)/2 and r t * = -/« 2 (p + 3P)/2. 



IV. SUMMARY AND DISCUSSION 

In this work we have derived the structure equations for hydrostatic equilibrium of spher- 
ically symmetric systems in a family of Palatini theories of gravity where the Lagrangian 
is a function of the form f{R ) R fJLV R^ lv ). Therefore, the TOV equations of GR have been 
extended to a much wider family of theories of gravity. We have shown that the differential 
character of the GR equations is retained, since we only have first-order derivatives of the 
functions ty(r),M(r), and P(r), but new nonlinear contributions (specially of pressure 
terms) appear due to the non-trivial role played by the matter sources in the determination 
of the connection. The corresponding limits to the cases of f(R) theories and GR have 
been explicitly computed and some typos on previous literature have been corrected. 



We note that with the expression for P rr obtained in Section flHI Ej) . one can further 
manipulate f 1 3 . 3 j) to arrange all the M r terms together on the left-hand side. That is the 
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form of the equations that should be used in numerical computations. However, since 
that representation does not provide any new physical insight, we have omitted that step 
here. We also note that the (unique) exterior solution of the structure equations is of 
the Schwarzschild-de Sitter type. Since outside of the star p and P vanish, we find that 
o~i = o"2 =constant, which implies that Q = S = o\ and T r r = Tt =constant, and the 
equations boil down to ip r = and M r = (f / 4af)\ vac r 2 , where (//4cr^)|„ ac is evaluated in 
vacuum and plays the role of an effective cosmological constant. 



Though interior solutions will be studied in detail elsewhere for different choices of the 
function f(R, R^R^), let us briefly discuss the potential effects of considering the following 
family of quadratic models: f(R,Q) = R + aR 2 /Rp + R^ V R^ V /Rp, for which R behaves 
exactly like in GR, R = —k 2 T, and Q = R ilv R! iV is given in (I2.16p . Here Rp is assumed to 
be some high curvature scale, such as the Planck scale. If we take a = —1/2, we find(23l. 25 



2k 2 ( P + P) 2k 4 { P -3P) 2 _ / An 2 (p + P) 
Rp 3Rp y Rp 



(4.1) 



At low energies, this expression recovers the GR limit, Q ~ (3P 2 + p 2 ) + 3 ^Rp + ^mt^ ~^ 
. . ., but at very high energies, positivity of the argument in the square root of (14. ip implies 
that K, 2 (p + P) < Rp/4, which clearly shows that the combination p + P is bounded 
from above regardless of the symmetries or particular configuration of the fluid involved. 
The interior solutions of this model must be very similar to those of GR except at the 
innermost regions of extremely compact objects, where the modified dynamics and the new 
pressure gradients should play an important role. Since the differential equations have the 
same degree as those of GR, we do not expect new solutions which may depend on free 
parameters, as it happens in other types of modified theories which introduce higher-order 
equations. Rather, the solutions of our set of equations must represent deformations of 
those found in GR. In fact, the GR solutions should be recovered smoothly and in a unique 
way in the limit Rp — > oo. The fact that this family of models leads to bounded density 
and pressure raises a natural question: can we find static solutions corresponding to objects 
denser than the black holes of GR? This question is pertinent because in GR there can 
not be static solutions with r — 2M(r) < 0, since they unavoidably lead to gravitational 
collapse and the divergence of energy density and curvature scalars. The fact that p and P 
are bounded in this theory suggests that such static solutions could exist. Exploring this 
possibility will be the subject of future research. 
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